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Abstract 

We study the the high spin expansion of the anomalous dimension for long operators belonging 
to the sl(2) sector of N = 4 SYM. Keeping the ratio j between the twist and the logarithm of 
the spin fixed, the anomalous dimensions expand as 7 = f(g,j)lns + f^°'(g,j) + 0(l/lns). 
This particular double scaling limit is efficiently described, up to the desired accuracy 0(ln 
in terms of linear integral equations. By using them, we are able to evaluate both at weak 
and strong coupling the subdeading scaling function f^(g,j) as series in j, up to the order 
j 5 . Thanks to these results, the possible extension of the liaison with the 0(6) nondinear 
sigma model may be tackled on a solid ground. 
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1 Introduction 



The set of composite operators 

Tr(V s Z L ) + . . . , (1.1) 

where D is a covariant derivative acting in all possible ways on the complex scalar field Z, consti- 
tutes the so-called si (2) sector of M = 4 SYM theory. The integer numbers s and L are called spin 
and twist, respectively. In the framework of the AdS-CFT correspondence [1] this set of operators 
received particular attention, also because of its connection [21 [3] to twist operators in QCD. In 
the high spin limit anomalous dimensions j(g, L, s) of (11.11) show a logarithmic divergence, 

7 (<7, L, s) = In s f(g) + f c {g, L) + o(s°) . (1.2) 

The function of the coupling f(g), which equals twice the cusp anomalous dimension of light- 
like Wilson loops, is known also as universal scaling function, where the term 'universal' takes 
into account its lack of dependance on the (fixed) twist L. This property, in particular, must 
imply that f(g) is free from wrapping effects, which makes the asymptotic Bethe Ansatz equations 
[H El El U\ predictive for its exact determination. On the other hand, the function f c (g, L) is not 
'universal', but, since it is related through functional relations [SI [9] to f(g), it contains a universal 
part. Finally, one has to mention that both / and f c are related [10] to form factors in scattering 
amplitudes. The function f(g) can be found from the solution of a linear integral equation [Til E] 
derived from asymptotic Bethe Ansatz. This equation was solved both at the weak [HI [7] and, more 
importantly, at the strong coupling [121 Ell [33] limit. Moreover, the function f c (g,L) was studied 
in [15] (weak coupling), [9] (strong coupling from string theory) and in the contemporaneous paper 
[16l[T7] (strong coupling, using linear integral equations derived from gauge theory Bethe Ansatz). 

The structure (11.2j) of the high spin anomalous dimension is preserved if one introduces the 
following scaling limit, 

L — 2 

s — > +oo , L — ► +oo , j = — fixed , (1.3) 

Ins 

proposed^ in [18] for the strong coupling (string theory) case and then in [15] for the weak coupling 
J\f = 4 SYM theory. Similar limits were also studied in one loop Af = 4 SYM - in [19] , when ln ^ L ^ 
is fixed - and in the strong coupling M = 4 SYM [61 [20], [2H [22] and string [191 [231 [21] cases, when 
, f L IT , is fixed. 

\n(s/ L) 

As we will show, in the case of limit ( 11.31) . relation ( 11. 2ft is replaced by 

oo 

1(9, J, s) = f(g,j) Ins + f^(g,j) + £ f^\g, j)(lns)" fc + 0((lns)-°°) , (1.4) 

k=l 



lr To be precise the limit considered in [TH [T5] is 



-oo , L — > +oo , j = - — fixed . 

In s 



Referring to expansion (|1 .4|) . the different limit (|1.3p docs not affect f(g,j), but gives easier forms for f^°\g,j). 
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where the notation 0((lns)~°°) stands for terms going to zero faster than any powers (lns)~ fc , 
k G N. Importantly, we may conjecture that wrapping effects can affect none of the functions 
appearing in (jl.4p . which thus are still completely determined by the asymptotic Bethe Ansatz 
equations or (non-)linear integral equation [TTj . One reason for that is the matching comparison 
with the string expansion of [9]. 

At small j the structure of the function f(g,j), called often generalised scaling function, 

oo 

f(g,j) = Y l fn(g)j", (i-5) 

n=0 

was investigated at weak coupling [J5] and at strong coupling [27J EE]- The generalised 

scaling function f(g,j) for j <C g was also shown [26] to coincide with the energy density [29] of 
the nonlinear 0(6) sigma model embedded into AdS§ x S 5 , thus confirming the related previous 
proposal [18] by Alday and Maldacena. Results for j ^> g are also available [30] . 

The aim of this paper is to investigate the structure of f^°\g,j) at very small j, both at the 
weak and the strong coupling limit. We expect f^°'{g,j) to be an analytical function of j, so that 
the expansion 

oo 
n=0 

will capture all its properties. Even if we report explicit results on functions fn\g) for n < 5, we 
will show that in general the various fn^ (g) can be obtained using a linear integral equation - giving 
the density of Bethe roots and of the so-called holes - which is equivalent to the asymptotic Bethe 
Ansatz equations in the high spin limit if one neglects terms of order (In s) -1 . This equation follows 
from previous results, mainly contained in [3TJ [32], and generalises the one we used previously 
[25j[27J[28] to determine the components f n (g) of the scaling function f(g,j) and which is equivalent 
to the 'FRS' equation, proposed in [T5j . 
The plan of this paper is as follows. 

In Section 2 we report on the one loop results. We show that the density of Bethe roots and 
holes can be found by solving a linear integral equation, which is exact if one neglects terms of 
order 0(1/ In s) and we give the expressions for fn (g = 0), when n — 1, ... ,5. 

In Section 3 we write the linear integral equation satisfied by the density of roots and holes at 
arbitrary values of the coupling constant. Using such density, we perform weak coupling compu- 
tations and give fn\g) up to three loops and up to n = 4. 

In Section 4 we re-write the linear integral equation for the density as a set of linear systems: 
from the solution to the n-th system one can get the n-th component in the expansion (11.61) . The 
solution to the n-th system depends - very similarly to the systems describing the components 
f n (g) of the generalised scaling function - on the previous ones. This reformulation of the problem 
simplifies the analysis of the strong coupling limit, which we do up to n = 5. 
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2 One loop results 

We begin our study by considering the sl{2) sector of Af = 4 SYM theory at one loop. In order to 
distinguish them from the exact (i.e. all loops) correspondents, all the one loop quantities (such 
as the density a and the roots- holes separator c) will be denoted by an index 0. As stated in the 
Introduction, we are interested in the limit (11. 3p 



-oo , L —> +oo , j 



L-2 
Ins 



fixed . 



In this limit, the density of Bethe roots and holes at one loop in perturbation theory, ctq(u), satisfies 
the linear equation 



a (u) = L 



ip ( iu J + ib I — V iu 



— ib 1 — iu — i— ■= 
\/2 



ib 1 — iu + i—= 
V2 



+ 



ib 1 + iu — i—= 
dv 



ib 1 + iu + i—= 
V2 



21n2 + 



— [^(l — iu + iv) + ib(l + iu + iv)} a (v) + O ( - — 



2n 



Ins 



(2.1) 



which, as we wrote, is exact if we neglect terms of order 0(1/ Ins). Further, we have to impose 
the condition 

1 



en 



duo~a(u) 



-2wj Ins + O 



Ins 



(2.2) 



which comes from the fact that the parameter cq defines the interval [— Cq, Cq] in which the L — 2 
internal holes concentrate, 

It is convenient to use the Fourier transform, ao{k), which satisfies the linear equation 



a (k) 



-2tt- 



f (1 



_]fcj 
e 2 



ifci 
e 2 



l-cos^| 



e 2 
sinh ^ 



sinh 

"°° dh 

vn 2tt 



\k\ 



AnS(k) In 2 



sin(/c — h)co sin hco 



k-h 



h 



O 



1 

In s 



(2.3) 



and the condition 



2 / — ff * 



sin fcc 



-27rj Ins + O - — 
Ins 



We can now give justifications of these formulae. We start from (3.52) of 



(2.4) 

where the L — 2 in 



L-2 



the last term of the first line is replaced by a sum on the internal holes Yl e lkUh . This equation, for 



h=i 



j^F (u) = <Jq(u) + 0(l/s), is exact if we neglect terms of order 0(l/s). Now, in order to express 
the sum on internal holes in terms of the density o"o(m), we use results from Appendix A. Referring 
to formula ( 1A.3I) . we realize that, if one neglects terms of order 0(1/ Ins), such sum is given by 
the first term in the right hand side of flA.3[) . which is linear in the density (Jq(u) = -^Z (u). After 
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doing this, we get expression (12.31) for the equation satisfied by the density of Bethe roots and holes 
in the Fourier space ^| . In a completely analogous way, when we express the counting of internal 
holes in terms of cr (u), we get (12.21 12.4p . 

As a further verification of the correctness of our starting equations, in the first part of Appendix 
B we prove that equation ( 12.11) is compatible with the corresponding equation - valid for very small 
u - that can be deduced from the results of [T9] . 

Having established on quite firm ground our starting point, we pass to analyze in detail the 
behaviour of the various quantities in the limit ( 11.31) . Consistency considerations coming from the 
analysis of ( 12.31 12.41) imply that in such limit the parameter Co, depending, via ( 12.31 12.41) . on Ins 
and j, expands, when j <C 1, as 



Co 



n=l v ' n=l x v ' ' 



Expanding also the condition (12. 4p we have 



-2irj Ins 



[2.5) 



(2.6) 



where we use the notation <7o(0), cro(O) to indicate the values at u = of the function cr and its 
second derivative, respectively. For such quantities expansions similar to ( 12.51) hold, e.g. 



OO OO / . 

a (0) = lnsg>^(0)j'1 + E4 M (0)i1 + 



n=0 



n=0 



[2.7) 



Explicitly we have 



or (0) = [-4 - 4j In 2 + 0(f)] In s - [8 In 2 + A lE + 0(f)] + O ( -L 

Ins 



and also that 



^'(0) = [56C(3)+ 0(j)] + 



1 

In s 



:2.8i 



(2.9) 



Now, inserting the expansion ( 12 .5[) in the condition (12.61) and using also (12.81) we can start finding, 
by equating equal powers in j and In s, the various coefficients Cg ' 
We get, without much ado 



5 c o 



,(o,i) 



7T 

4 

71 



.(0,2) 



7T 

4 



In 2, 



.(0,3) 



7T 



(ln2) 2 



.(i,i) 



(21n2 + 7£ ) 



C ( 1 < 2 ) = ^ln2(21n2 + 7i3 ), 



|7r(ln2) 2 (21n2 + 7£ ) + I ^7r 3 C(3), 



in such a way that one can write 



Co 



-j --\n2f 

4_ J 4 J 



7T 



7T 



+ '- In 2 (2 In 2 + 7s ) j 2 



(21n2 + 7B )j 
7 



i7 r(ln2) 2 (21n2 + 7j5 ) - ^^(3) ) j + 0(j 



T— + 

ms 



2 Actually, in order to get equation (|2.3|) . we have to use once identity (|2.4 



(2.10) 
(2.11) 



(2.12) 



(Ins) 
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All this information can be used to compute the anomalous dimension at one loop 7^2 up to the 
desired order in j. From the formulas - exact if we neglect terms O (^7) - 



du 



7 9 2 =gE 0l E = - —e(u)a (u)+ —e(u)a (u 



dk 
4^ 



e(k)a (k) + 



2tt 

dk 
4^ 



du 



CO 



2vr 



e(k) 



dh 
2^ 



00(h) 



^sin(/c — h)co ^smhcQ 



k-h 



h 



(2.13) 
- lnsje(O) , 



where E denotes the energy of the spin —1/2 Heisenberg chain and the function 



e(u) 



u 2 + 



e(k) = 27re 2 ; 



(2.14) 



we get, expanding in powers of j: 
= g 2 In s 



+ 9 1 



4 7j 



4-41n2j + 
7C(3)vr 2 



7C(3)vr 2 - 



21 



7C(3)vr 2 ln2 . 4 
12 3 



77r 2 (ln2) 2 C(3) 31tt 4 C(5) 



640 



f + 0(f 



12 



(2 In 2 + lE ) f + 7C(3) f ln2 (2 In 2 + lE )f 



77r 2 (ln2) 2 C(3) 31tt 4 C(5) 



160 



4 

(2 In 2 + 7s ) j 



960 



As we show in Appendix B, one can also get such an expression improving the results of [19] 
in order to include also contributions of order (lns)°j fc (as far as we know, authors of [19] are 
interested in order lnsj fc contributions and develop their calculations accordingly). This reinforce 
our belief in the goodness of ( 12. II) as description of the leading and subleading high spin behaviour 
in one loop si (2) sector of M = 4 SYM. 



3 All loops equations 



Let us now pass to study the anomalous dimension as a function of the coupling constant g. As 
usual, we split the density of roots and holes, a(u), as a{u) = a (u) + a H (u), where a (u) is the 
one loop contribution and (Th(u) is the higher than one loop contribution to the actual density. It 
is convenient to introduce the quantity S(k), related to the Fourier transforms of cr (n) and (Th{u) 
as 



, N 2sinhf 
S(k) = ^ 



a H {k) - 



- 1M 
e 2 

sinh 



+ °° dp A . sin(A; — p)co e '2 
<ro(P) ; + 



2tt 



k — p 



sinh 



+ °° dp * / \ sin(A; — p)c 
a(p) 



2tt 



k — p 



(3.' 



where c indicates the separator between the holes - concentrated in the interval [— c, c] - and the 
Bethe roots. 
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The function S(k) satisfies the linear equation 



S(k) = -^[l-J {V2gk)} 



7T\k\ 



+°° dh 

oo \h\ 



^r{-l) r+l J r (V2gk)J r (V2gh)- 



sgn(kh) 



+ 



oo oo 



r=2 u=0 



+ sgn(/i) ^^c r , r+l+ 2u{-l) r+u e-^ (j r ^(V2gk)J r+2u (V2gh) - J r ^{y/2gh)J r+2l ,{y/2gk) 



JAl 
e 2 



{ ^Ls(h) - An In 2 - ttj In s^— ^ 
•> sin h J — l sinh - — ^ 

sin(/i — p)c sinpc 



2vr- 



2 COS t5§ 



sinh 

1M 
e 2 



sinh ^ 



sinh 



Z7l 



h 



V 



p 



Ins 



(3.2) 



And this relation has to be solved together with the condition 

' +oc dk sin kc 



-2nj In s + O [ - — 



1 

In s 



(3.3) 



2tt~ v "' k 

To justify equation (I3.2p we start from (4.11) of [32], where in the last term of the first line L — 2 



L-2 



is replaced by a sum on the internal holes e an d the in the arguments by u%. Moreover, in 

h=l 

the integral in the second line we can replace the extremes ±6q with ±oo and -^F (v) with <Jq(v) 



( 12. II) . The equation obtained in such a way for -^F H (v) = cth(v) + 0(1/s) is exact when neglecting 
terms of order 0(1/ s). Now, we can use results in Appendix A §, which are valid for a general 
counting function Z(u) such that Z(c) = —ir(L — 2) + 0(1/ Ins), which is exactly condition (13. 3p . 
Using (1A.3|) . we evaluate the sum over internal holes keeping only the first (linear) term in the 
right hand side of this formula, since we want to neglect 0(1/ Ins) contributions. After passing to 
Fourier transforms and defining ( 13. ip . we get eventually equation (13.21) . 

Finally, we notice that the relation of the function S(k) with the anomalous dimension 7(g) is 
a generalisation of an identity found in |33j : 

1(9) = 25(0) . (3.4) 

The proof of such identity follows by the simple comparison between (13.21) and the expression for 
the anomalous dimension 



1(9) 



du 
2^ 



q 2 (u)a(u) + 



du 
27 



q 2 (u)a(u) , q 2 (u) 



x + (u) x~(u) 



X 



u) 



u± 1/2 




(3.5) 



(u± 1/2) 2 

A little care is needed in order to interpret (13.40 : the term proportional to g 2n , which depends on 
the density at n + 1 loops, gives the n-th loop contribution to the anomalous dimension. 

In the next subsection, we will compute at two and three loops the anomalous dimension in 
the limit (ll.3p . New results will be the subleading, 0(lns°) contribution, which we will evaluate 
as power series in j, up to j 4 . 



3 We remember that <j(u) 



■j^Z(u), where Z(u) is the counting function. 
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3.1 Two loops anomalous dimension 

Let us expand the linear equation for S(k) at order g 4 , i.e. at three loops approximation for the 
density. Before doing this, it is convenient to write the expression for S(k) at the order g 2 . From 
the general expression we get 



+00 



dh \[2gh \[2gk 1 — sgn(kh) _w 
\h\~^2 2 2 6 



7i j In s 



l-e¥ l-e-Vcos^§ 
-2tt 



e 2 



sinh ^ 



sinh ^ 



sinh *y 



dp 
27 



4tt1ii2 5(7i) - (3.6) 
sin ( h — p) Co sin pco 



P 



P 



Comparing (13 .4p with (13. 6p we get 



9- 



L \ k \ + o7q 



(3.7) 



Let us pass now to the three loops (order g 4 ) density. In this case the relevant equation reads 
-,4 



SAk) 



16 7T \k\ 



+ °° dh V2gh \/2gk 1 - sgn(kh) _hj n\h\ 

00 



e 2 



H ,+00 



sinh 



dp,, s 



2 2 2 
sin(/i — p)c sinpc 



sinh 



h — p 



P 



+ 



7T 



2 

n \h\ 



--g 4 k 2 h 2 



(V2gh) 3 V2gk (V2gk) 3 V2ghi 1 - sgn(kh) _w 

e 2 



16 



irj In s 



, M 
1 — e 2 

sinh ^ 



2tt 



16 



1 — e 2 cos ^| 
sinh ^ 



-4vrln2 



(3.* 



e 2 



sinh -rr 



<ip 
27 



o-o (p) 



sin(/i — p)co sinpco 



P 



P 



It follows that the two loops contribution to the anomalous dimension, 7^4, is given by 



7o* 



3_ 

2tt 



(i/ze 2 



IMf tt|^| 



1M 
e 2 



^ sinh 1M 9 v 7 1M 



dp 
27 



<r(p) 



sin(/i — p)c sinpc 



P 



P 



+ 



f- / dhh 2 e~^ 
8tt 



-Att In 2 5(^i) — 7Tj In s 



e 2 



sinh 



l'<l 



(3.9) 



2tt 



1 -JM hs 
1 - e 2 cos ^| 

sinh 



e 2 



sinh y 



+00 



dp 
27 



0b(p) 



sm(h — p)co sinpco 



h — p 



P 



It is not difficult to compute, neglecting orders 1/s, the contributions to 7^4 coming from the first 
term in the first line, the terms in the second line and the first term in the third line. For these 
contribution we get 



7 sM = -2s 4 C(3)(jln S + 2) - 9 -n 2 lg 2 + 6^ 4 C(3)j Ins - 2^(3) 



(3.10) 



-6^C(3)+4A(3)jln S -^7r 2 7 9 2- 



The computation of the contribution, which we indicate with 7^ 2, coming from the last term in 
(13. 9p relies on the one loop results. If one wants to restrict to terms containing powers of j not 
higher than 4, the relevant integral to compute is simple: 



g 4 f +OQ h 4 g 4 

^From the one loop results we have, up to the desired order 
cfco(O) = - 

r3 



62C(5)c^ (0) 



(3.11) 



16' 



:J 3 + — J 4 In 2 



In s 



+ 



7T 



7T 



(2 In 2 + 7_b) j — 3— j In2(21n2 + 



+ 



In s 



(3.12) 
(3.13) 



in such a way that the corresponding contribution to the energy reads 
,4 



7 9 4 ,2 



9_ 



62C(5) 



16 



3 ^3 

f + — f In 2 



Ins 



TV 



.3 3 

62C(5) — (2 In 2 + lE )f - 3^/ In 2(2 In 2 + lE ) 



8 



O 



In s 



(3.14) 



We are now left with the calculation of the contribution - 7^ 3 - to the anomalous dimension 
coming from the last term in the first line of (13.91) . Since we restrict to terms containing powers of 
j not higher that 4, at desired order such contribution equals 



%*3 



a 2 r°° h 2 14 

-Ufa®]* / dh— tt = -/-C(3)[cV(0)] 92 . (3.15) 
D7T a J smh I 37r 

We have now to compute <r(0) and c up to the order g 2 . For what concerns cr(0) we have 

a (0)^[-A-Aj\n2]\ns-[8ln2 + A lE } + [a H (0)} g 2 + o(^-J , (3.16) 
where, from (13.11 13. 7[) . one has 



[Mo)L 



7T 



2 sinh 



l'<l 



[S(/i)], 2 + 0(j J ) = 14(j In a + 2)^(3) + ^ +0[ — \. (3.17) 



In s 



On the other hand, the condition (13.31) for c at the desired order reads as 

2a(0)c = -27rjlns . 
Using such relation, we can get c at the order g 2 and up to the order j 2 : 



(3.18) 



7T 



16 In s 



j(l-21n2j) 



7T 



( 7£ + 21n2)(-2j + 61n2j s 



16(lns) 2 

[28£ 2 C(3) + y 7g2 + 14j(? 2 C(3) Ins] + O (J-} . 



(3.19) 



Coming back to the energy, we have 



14 

V,3 = (? 2 — C(3)vrjln S [c 2 ] 5 2 +0(j 5 ,lns) 
which, after using the one loop expression for cq up to the order j 2 , 



Cq 



71 



7i 



In 2/ • Oif 



+ 



7T 



- In 2(2 In 2 + lE ) f + 0(f) 



(21n2 + 7fi )j + 
1 



In s 



O 



(Ins) 



(3.20) 



(3.21) 



gives 



7^,3 



^C(3)AV[(1 - 3j In 2) + glgj±3g (-3 + 12 In 2 j) 



In s 



7T 



[2% 2 C(3) + yV + 14j^ 2 C(3) Ins] + O 
Summing together (13.101 13.141 I3.22p we get the energy at the order g 4 



(3.22) 



lg 4 



VC(3) + 4A(3)jln S -|-7r 2 7 fl2 + 



— 62C(5) 

7T 



7T 

16 



7T 



f + —f In 2 



In s + 



62C(5)(21n2 + 7£ ) 



7T 

+ ^C(3) g 2 n 2 f [(1 - 3j ln2) + 2 ^^ (-3 + 12 j In 2) 

2 / 1 

Y(31 4- — 'v.a 4- 147 a 2 f(?>\ In si + O I 

In s 



[28(? 2 C(3) + y V + WC(3) Ins] + O ( — 



Alternatively, we have 
r 2 

lg 4 = 

31 



7T 

"6 



— s' 2 v + g 4 m s 



4C(3)j-yvr 2 C(5)j 3 + lA(3)j 3 + 



49 7 
+ y vr 2 In 2 C(5) j 4 + -tt 2 C(3) 2 j 4 - -vr 4 C(3) In 2 j 4 + 0(f) 



+ 



r 31 

^ 4 [-6C(3) + T 7r 2 C(5)(21n2 + 7i; )f 

7 a , „ „ x ./i 93 



49 
12 



2 -3 



7T 2 C(3) 2 J 



7 

12 



(3.23) 



A(3)(31n2 + 7s )j ; 



+ -tt 4 C(3) In 2(5 In 2 + 2 lE ) f - — tt 2 C(5) In 2(2 In 2 + lE ) f - 
4Q 

- -^C(3) 2 vr 2 (41n2 + 7£ )j 4 + 0(/ 



In s 



(3.24) 



3.2 Three loops anomalous dimension 

For the three loops anomalous dimension, the calculation follows the same route as in the two 
loops case. We simply report the final result. The three loops anomalous dimension 7s e expands 
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as 

oo oo 

7g6 = In S ^ fnj f + fnl* I" + °(V 111 S ) ( 3 - 25 ) 
n=0 n=0 

where the coefficients f n>g e were already computed [15] (at least, explicitly up to n = 4) and the 
new coefficients 8) still up to n = 4, read as 

= 20C(5) + |tt 2 C(3) + l X^ E (3.26) 
/S = /J = (3.27) 
/5 = -^7r 6 C(3)(21n2 + 7E ) + ^A(3)(421n2 + 49 7E ) + |7r 4 C(5)(21n2 + 7B )- 

- yA E ((5) + ^ AO) 2 - In 2 + 7 s) - ^ 2 C(3)C(5) (3.28) 
= ™7r 2 ln2(21n2 + 7£ )C(7) + ^ 7 r 2 (21n2 + 7£ )C(3)C(5) + ^vr 2 ln2C(3)C(5)- 

- ^tt 4 In 2(2 In 2 + 7s )C(5) - ^vr 4 (ln 2) 2 ((5) - ^tt 4 (2 In 2 + 7£ )C(3) 2 - ^tt 4 In 2((3) 2 + 
+ ^ 2 C(3) 3 + ^tt 6 In 2(2 In 2 + 7s ) C (3) + ^ 6 (ln 2) 2 C (3) . (3.29) 

4 Systematics of the subleading term 

We now want to put the linear integral equation (I3.2p in a form which is more suitable for analysis 
of both the weak and the strong coupling limit. In brief, instead of working with one linear integral 
equation (I3.2p . we will be concerned with (two) linear infinite (i.e. containing infinite equations) 
systems. 

As far as the dependence on Ins is concerned, equation (13.21) and its solution split in a part 
proportional to Ins and a part proportional to (Ins) . 

We concentrate on the latter, which we call S^°\k), since the former has been extensively 
studied pH[28]. 

We now restrict to the domain k > and expand S^°\k) in series of Bessel functions 



S(°\k) = f2sl°\g) Jr{V ? 9k) . (4.1) 



r=l 



As a consequence of (13.21) . the coefficients Si -(g) satisfy the following system of equations 

dh J 2p ~i{V2gh) 



SUM = 2^7^,1 + 4(2^-1) f f J ^ ( I f *° + AfUa) - 2(2p - i) E Zt-iMSgig) 

771=1 

(4.2) 

Si°j(9) = 4 + 8p r f J2 ^f ] + + Ap Z^.^S^g) - 4p Z^(g)S®(g) . 

777=1 771=1 
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where, as usual, we introduced the notation 



dhJ n (V2gh)J m (V2gh) 

Z n ,m(g)= I y • ( 43 ) 







and where we introduced the function J2 P -i{x) which coincides with the Bessel function J2 P ~i{x) 



for p > 2 and with J\(x) — | when p — 1. In (14. 2 p the terms Ai°\g) is the term proportional to 
(Ins) of the quantity: 

sin(/z — p)c sinpc 



Jo 27i7i sinh| J.^ 2tt 



/i — p p 



Hp) ■ (4.4) 



Eventually, anomalous dimension at order (Ins) , 7 < -°\ is extracted from S[°\g) by means of the 
formula 

7 (0) = V2gSf\g). (4.5) 

To be complete, we write down also the equations satisfied by the part of S(k), which is linear in 
Ins. Let us call such part S^^k). Expanding in series of Bessel functions 



S^\k) = f2s^\g)^^, 



k 



(4.6) 



we find the following system of equations for the coefficients 

S£>M = 2V2 g d p ,i ~ 2(2p - l)j f°° f ^ + A£\( 9 ) 

Jo n ea + 1 

oo 

- 2{2p-\)Y,Z2 P -iAg)s { n ; 1 \g) 



m=l 



(4.7) 



oo 



Jo h> + 1 

oo oo 

+ Ap Z 2p , 2m ^ 1 (g)St 1 l 1 (g) - Ap £ 

m=l m=l 

where A^f (g) is the term proportional to Ins of the quantity (I4.4p . 

This system has been used in the series of papers [251 123 [28] in order to compute the (strong 
coupling limit of) the generalised scaling function f(g,j). In the next subsections we will adapt 
the steps of [25], [271 EE] to system ( 14.21) in order to study the subleading function p°'(g,j) as series 
in j, up to the order j 5 . 

4.1 Slicing in powers of j 

Since we are in the limit (11. 3p . the function S^°'(k) admits an expansion in powers of j, 

oo 

S {0 \k) =Y,S (0 ' n \k)j n . (4.8) 

n=0 
12 



Consequently, this way of expanding extends to the coefficients Si°\g) (14.11) as well: 

oo 

Si°\g) = Y,Si > n) (g)f. (4.9) 



n=0 



In expansion (14 .9p the coefficients Si°'°\g) - which give the part of the function S^°\k) inde- 
pendent of j - satisfy the system of equations 



Si°J°\(g) = 2V2g lE 5 p , 1 + 4(2p-r 



dh J 2p -i(V2gh) 
~h e h - 1 



(4.10) 



m=l 



S£' 0) (S) = 4 + 8p 



dh J 2p (y/2gh) 



h e h 



m=l 



m=l 



This system of equations is studied in the contemporaneous paper [T7] - the coefficients Sr°'°\g) 
being there denoted as S e r xtra (g). Therefore, we pass to study the function Sj°' n \g), when n > 1. 
For simplicity's sake we limit to the cases n = 3,4, 5 (for n = 1,2, Sr (g) = 0). The relevant 
forcing terms are obtained after developing (14.41) for small c up to c 5 : 



A^ n \g) 



+oo 



dh J r {V2gh) 
2nh sinh | 



-\h 2 c 3 a(0) + ^ 4 cV(0) - -^ 2 cV 2 (0) 
3 60 10 



(Ins) 



Now, from the expansions in powers of In s: 

cr(0) = Ins a { ^ 1} (0) + a (0) (0) + O f — 1 

\ In s J 

we get (3 < n < 5): 

" +oc cZ/i J r (V2gh) r 1 



3 < ra < 5. 
(4.11) 



cC) + (In a)" 1 cW + O^^j) , (4.12) 



4°'%) 



27r/i sinh | 



^ c (0)V°)(0) + 3c(°) W^O) + 



(4.13) 



60 



+ JLtf ( C (0)V(°)(0) + 5c<°> W^O)) - -U 2 fc(°) 5 4 0) (0) + Sc^V 1 ^ 1 ^) 
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.(o) 4 Ji)J-i), 



J 3 



Keeping in mind the expansions in powers of j: 

oo oo 

c (fc) =^c M j n , a (fc - 1} (0) = J2 aik ~ 1,n) (°)j n > fc = 0,l 



(4.14) 



n=l 



n=0 



we can now specialise (14.131) to the cases n = 3, 



4°' 3) (s) 



dh J r (V2gh) 
2nh sinh | 



/i 2 (c^V^O) + 3c (0 ^ 2 c^^- 1 '°\0) 



(4.15) 
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n = 4, 

A^\g) = r 



+00 



dh J r (V2gh) 
2nh sinh I 



_i/,2 3c (0,l)^(0,2) a (0,0) (0) + 3c (0,l)yi,l) a M,l) (0) + 

3 V 



+ 3c^ 2 c^a^°\0) + Qc^c^c^a^°\0)) (4.16) 



and n = 5: 

4°' 5) G?) 



+ °° ^_ J r(V2gh) r_ W 3 c (o,i) 2 c (o,3 )(J (o,o) (0) 3c (o, 2) 2 c (o,i )(T (o,o) (0) 
27r/i sinh I L 3 V v 7 v 7 

+ 3c( ' 1 ) 2 C ( 1 ' 3 V^ 1 '°)(0) + Gc^c^c^V-^O) + Bc^'V ' 3 ^ 1 '^*- 1 ' ^) + 
+ 3c^ 2 c^a^°\0) + Sc^'c^V" 1 ' 1 ^) + 6 C ( ' V ' 2 ^ 1 ' 1 ^^ 1 ^)) + 

+ l/ l ^ c (o,D 5 a (o.o) (0) + 5 C (°' 1 ) 4 C ( 1 ' 1 ) ( t(- 1 ' )(0) N ) - 
60 V / 

_ L h * (cm 5 a m {0) + 5 C (o,i) 4 c (^)4-^(o)) 



(4.17) 

On the other hand, the c can be expressed in terms of the <r(0) by means of the condition (13. 3p . 
At the relevant order in c such condition reads 



2a(0)c+ -ct 2 (0)c 3 = -27rjlns 
3 



The order In s° of such equation gives 



2(7 (-i)(0) c ( 1 ) + 2a< >(0)c<°) + ^(O^V 1 ) + -4 0] (0) C W 3 = 



Specialising such equation at order j, we get the condition 



(o,o) (o) 



(T (-i,o)( ) 



.(0,1) 



a 



7T- 



(0,0) (o) 



(J (-1,0)(0)]2 • 



At order j 2 we have 



.(1,2) 



a 



(-1 



D(0) 



,(i,i) 



<t(°- >(0) 

Going at order j 3 we get 

a^fm (1)2) _ ^ ( °' 0) (0) 
o-(-i,o)(o) C ^1,0(0) 
^(o.o)^ [(T (-i,i)( )]2 

3 



c (0,2) = _ 2n 



a (o,o) (Q)(T (-i,i) (0) 

^(-^(o)] 3 



(4.18) 



(4.19) 



(4.21) 



(4.22) 



.(1,3) 



1 rr M ' 0) 
(0,3) _ til 



(0) 



1 rr (0 ' 0) fm 



2 ( j(- 1 '°)(0) 



,(0,1)^(1,1) _ 



.(o,i) c 



6 cr( _1 '°) (0) 



3tt- 



2 7r 3^ (0 ' 0) (0)^- 1 '° ) (0) , tt 3 



t7(-i,o)(o)]£ 



6 [a(- 1 '°)(0)] z 



4 We remember the formulae (relations (5.18-5.20) of [25]): 



„(o,i) _ _. 



tr(-l.°)(0) 



[ a (-l,0)( )]2 ' 



.(0,3) 



.3 „("1,0) 



(0) \o 

— 7T 



(-1 



1} (0)] : 



6 [<r(- 1 ,°)(0)] 4 [<t( _1 ,°) (0)] 3 



(4.23) 



(4.20) 
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In such a way, we obtain 

A^\g)=r 

and, also, 



+00 



dh Jr{V2gh) ( 2 
2i\h sinh | 



-tt 3 /* 2 



g(°'°>(0) 

(-1,0)(0)]3 



dh J r (V2gh) 3,2^ (0 ' 0) (0)^ { ^ 1 ' 1) (0) 
Lix h — 



2nh sinh 



(-i,o)(o)]^ 



(4.24) 



(4.25) 



For what concerns the term proportional to j 5 , we get, after some calculation: 



4°' 5) G?) 



+00 



dh J r (V2gh) 
2nh sinh | 



7T 



r (o,o) 



(0) 



5 [trf" 1 ' )^)] 



+ 12tt 



, a (o,o )(0)[(T (-i,i)(o)]2 (o,o) (0)(J (-i,o) (0) 

' r 7 1 , 7T7 VT r / . ^ ; 777 



(-1,0) 



(0)] E 



(-1.0)(0)]6 



I 1,4 5 ^ 0) (0) 

15 [^-^(O)] 5 ] 



(4.26) 



Now, in analogy to what done in [28], we introduce the 'reduced' coefficients Si k \g) which satisfy 
the system (4.23) of [28], i.e. the 'usual' system with the BES kernel and with forcing terms 



+ °° dh h2k ^ J r (V2gh) 
2ir sinh | 

(0) 



(4.27) 



Using notations of [28], we can write the various fn', n = 3,4,5, in terms of S[ k \g) and of the 
density and its derivatives in zero 

2 _ 



fl%) 



- 



[ a M,o)( )]4 

.(0 
2 



(4.28) 
(4.29) 



1 /tt 5 4°' 0) (0) 



3 V 5 
7T 5 a(°'°)(0) 



12tt : 



ig (0,0) (Q)[a (-l,l) 

[a(-W(0)]t 



— IT 



i( T( -°)(0)^- 1 - 0) (0) 

[ (T (-1,0)( )]6 



V2^; i) (^) + 



(4.30) 



15 [^-^(O)] 8 

For comparison, corresponding formulae for the coefficients f n (g) of the function f(g,j) are, (5.30- 
32) of [28], 

7T 3 1 



h{g) 
h{g) 



3 [a(^.°)(0)] 2 

r(-l.l)(0) 



2 3 (7 

"3^ [^=W(0)p 



(4.31) 
(4.32) 



(-14) 



7T 



7T a, 



(-1,0) 



(0) 



V2g~S^\g) 



7l 



V2gS[ 2) (g).(4.33) 



60 [at-^HO)]' 



[ a (-l,0)( )]4 1 5 [c7(-l,0)(0)]5. 

Expressions (I4.28[ I4.29[ I4.30P interpolate from weak to strong coupling and for the moment lack of 
an explicit form as functions of g. Their weak coupling expansion in powers of g 2 were given, up 
to g 6 , in Section 3. In the next subsection, we will study and explicitly find their strong coupling 
limit. 
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4.2 Strong coupling 

It is interesting to consider the strong coupling limit of the equations obtained in the last subsection. 
For that purpose, we need to know the values of the density and its derivatives in zero and the 
quantities S[ k \g) as well. Results concerning a^~ 1,n \ i.e. the part of the density proportional 
to Ins, and S[ k \g) are reported in [28]. For what concerns the contribution to the anomalous 
dimension proportional to Ins, i.e. f(g,j), we remember [TBI 126] that in the limit (11.31) and when 
j <C g it coincides with the energy density of the 0(6) sigma model. The nonperturbative (infrared) 
regime of the sigma model, j <C m(g), where 

m{g) = kg*e~% + ..., k = ^Tzfy > ( 4 - 34 ) 

where the dots stand for subleading corrections, makes contact with the double limit j < 1, 
g — > oo, considered in this subsection (and also in [28]). In such regime the scale of the mass is 
given by m(g) and expansion of f(g,j) for small j and large g was successfully checked [271 128] 
against analogous non perturbative expansions [29] in the sigma model. 

It remains an open question if, as well as f(g, j), also (g, j) shows connections with quantities 
of the 0(6) sigma model. This is a further motivation which leads to study the strong coupling 
limit of f^(g,j), which however is an important problem in itself, since its results can be checked 
against string theory data. We concentrate on fn (g), for n = 3, 4, 5: as follows from (14.281 14.291 
14.301) , the calculation of their strong coupling limit can be finalised if we know the large g behaviour 
of the density and its second derivative in zero, o"( 0,0 ^(0) and a^'°\o). From results of [T7], we 
know that at large g and at the leading order 

^°)(^ = -ln^^- 1 '°)(^ + ... , (4.35) 

where Sp\g) is a solution of (14.101) and Sp~ (<7) of (14. 7p with j = (i.e. the BES system). 
From this relation we can deduce that, at the leading order, the density in zero reads as 

a JJ' 0) (0) = -\nga { - 1 ' 0) (0) + ... = - lng[A - Trm(g)] + . . . a m (0) = -4\ng + ... (4.36) 

Analogously, for what concerns the second derivative, we get 

4°' 0) (0) = 56C(3) + ... (4.37) 

since the higher loops contributions, ^ lng m(g), are exponentially depressed. 

Using the results contained in [28], we get the strong coupling relations (which are valid only 
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at the leading order), 



fl%) 
4%) 



h(g) 
h{g) 



81n# 
Txm(g) 
\2\ng 
Txm(g) 
16 lng 
ir 2 m(g) 1 

21ng 

Tim{gY 



it hig 
3m(g) 2 
\ng 



+ 



In 2 



ln2 + - 



m{gf 
2 7T 2 In g 



3m(g) 5 
16 Ing 



irm(g) 



V2gS[ l \g) 

h(g) 



In 2 



7r 3 In g 
30m(#) 4 



(4.38) 
(4.39) 



(4.40) 



We see that there is a simple proportionality between fn (g) and (n — l)f n (g) for n — 3, 4, which, 
however, is lost when n = 5. This could suggest a relation between f^°\g,j) and f(g,j), but at 
the moment our data are not conclusive about its form. 



5 Conclusions 

A systematic procedure for the computation of the subleading correction f^°'(g,j) in the high spin 
limit (11.31) to the anomalous dimensions in the si (2) sector of M = 4 SYM is developed. The 
method is analogous to the one used in [28] to study the leading term f(g,j) and is based on a 
linear integral equation describing the behaviour of the density and the observables at high spin. 
We first found one loop (subleading) contributions to the anomalous dimension (Section 2), finding 
agreement with results obtained by another method, which is an improvement of the techniques of 
[T9"] (Appendix B). Then, we performed weak coupling expansions (two and three loops, Section 

3) . Finally, we wrote the linear integral equation for the density as linear infinite systems (Section 

4) . This was particularly convenient, for it simplifies the study of the strong coupling limit of 
f(°\g,j). In this respect, we explicitly found (Subsection 4.2) the large g limit of fn\g), for 
n = 3,4,5, comparing their expression with f n (g). Although the string action may reduce to the 
0(6) non-linear sigma model, the question of its appearance in the sub-leading scaling function 
still stays as an open question. 
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A Estimate of the holes contribution 



We now give an estimate of the sum over the internal holes Uh of a generic function 0{uh). Results 
we present here are valid both at one and all loops cases. We start from the exact hole expression 
depending on the so-called counting function Z(u) §, 



L-2 



h=l 



du 



2tt 



0(u)Z'(u) + Im 



C ± (v-i6)^\n[l + 5e iZ ^} + 
„ n dv 



(a.i; 



+ Im 



^-0(-c + iy)4- ln[l + 5 e iZ ^ c+ ^\ + Im 
7r dy 



° dy O(c + iy)^-ln[l + 5e iZ( - c+i rt] 
dy 



7T 



The right hand side of (lA.lj) does not depend on e as far as no poles of the integrands lie in the 
region Ymv < e, |Reu| < c. The constant 5 is equal to ±1, depending on the parity of L. Keeping 
e small, but finite, we suppose that 

|e iZ(2) |<l, (A.2) 

when z belongs to the integration contour of (IA.1I) . In the large s limit this is justified, since Z'(v) 
is proportional to — In s. Using such approximation, we can replace all the ln[l + e %z ^] with e tZ ^ z \ 
Then, since the expression we get is still independent of e, we find convenient to evaluate it when 
e = @. We obtain 



g 0K) = _£* 0(uW « )+ , to £|o(.)^" 



(A.3) 



The nonlinear term in this expression 



m f C ± (v)^-e lZ ^ 

TT dv 



(A.4) 



can be estimated after the change of variable x = Z(v) (we remind that Z'iv) < 0) and after using 
the formula 



dxf(x)e° 



Em: 



j(k ){x) 



(A.5) 



fc=0 



where f^ h \x) denotes the /c-th derivative of f(x). We get 



NL = 51m 



7T * — ' Ox 

k=0 



Z(c) 



Z(-c) 



(A.6) 



The first terms of such series are 



NL = 5-0{c) sinZ(c) + 5-%^- cos Z(c) + 

7T 7T Z (C) 



(A.7) 



5 The counting function is connected to a(u) by the relation Z'(u) = <j(u). 

6 We can put directly e = since, after the approximation (|A.2[) . the integrand is regular on the real axis. 
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where the dots represent terms containing higher powers of Z'(c) in the denominator. Now, if c is 
chosen such that Z(c) = — 7r(L — 2) + 0(1/ In s) - which is exactly condition (12.41) for the one loop 
and (I3.3P for the all loops case, respectively - the first term in (1A.7I) is 0(1/ Ins). The second term 
is 0(1/ Ins) as well, since Z'(c) - for generic g - is proportional to Ins. After carefully evaluating 
all the terms in the series, one can show by similar reasonings that the k-th term is 0(l/(lns) fe ): 
this ensures the convergence of the series in (IA.6I) . which gives the nonlinear contribution to flA.3j) . 
Evidently, such nonlinear contribution is 0(1/ Ins). 



B Contact with 03 



We now show that our one loop results of Section 1 can be obtained starting from the equations 
of [19] . In this paper, authors are interested to the contribution to the energy proportional to In s 
(i.e. the part lns/(j), where f(j) is the one loop contribution to the cusp anomalous dimension). 
However, their equations can be used to compute also the subleading part, proportional to Ins , 
both at the level of equation for the density and at the level of explicit computations of energy 
eigenvalues. 

Let us start from the equations describing the distribution of roots and holes. In paper |19j . 
authors start from the equation 



25 n In s In 



L-l 



2 r(± + i5 n ) 



2^ r i 



3=2 



r (1 + i5 n — i5j 

iS n + iSj 



~k r , 



(B.1) 



which describe the distribution of internal holes, denoted as 8 n , 2 < n < L — 1. This equation 
suggests the definition of the counting function 



L-l 



Z (u) 



, T TU-zu) >v ^ T(l + iu-i6j, 
-4 In s u + iL In — )f '- + % > In — . 

^ r(i - in - is*) 

j=2 V 3 J 



rg + i«) 



(B.2) 



Indeed we have that Zo(S n ) = n(2n — L — 1) 



\L+1 



In terms of the derivative cr (u) = -^Z (u), we have the equation 



a {u) 



-4 In s + L 



( 2 ~ iu ) +^ \ \ + iu 



L-l 



^[-0(1 +iu- i5j) + ip(l -iu- iSj)] , (B.3) 



3=2 



which after expressing the sum on the holes in terms of the density of holes — ^o"o(w) turns out 
into 



&o{u) — — 41ns + L 



1 



en 



CO 



dv 
2^ 



ip [ - — iu)+if)[ — + iu] + I — ao(v)[i[)(l + iu — iv) + — iu — iv)\ . 

(B.4) 

Such equation describes the distribution of holes and, consequently, it is valid only for — cq < u < 
Co- And, indeed, one can verify that our starting equation (12. ip reduces exactly to flB.4p when 

—Cq < u < Cq. 
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After having shown that our equation for the density of roots and holes agrees with the cor- 
responding relation that can be deduced from results of [19], we pass to compare results for the 
eigenvalues of the energy: we verify that expression (I2.15P can be obtained starting from formula? 

of p]. 

Referring still to equation (1B.1I) . describing the distribution of internal holes, we develop such 
equation up to the order S^. We obtain the relation 



25 n \xxs-L^ I 1 -) 6 n+ ^ (l) ^ML-2)^l)5 n -^^\l)5l-\^\l)5 n ]T 6^0(5, 



i=2 



(B.5) 



which we rewrite as 



2 



2\ns — Lip (I) + (!) 5 2 + {L _ 2m) _ ^m{i)8i - W 2) d) <? + o{8% 

(B.6) 

Therefore, we can write that 



2 r 



+ o 



21ns - Lij) (|) +(L- 2)^(1) 
Si 



¥ {2) (I) *l - WW - l^ 2 (i) Etf + o(* 



21ns - Lifj (|) +(L- 2)^(1) 



(Ins) 2 



(B.7) 



^From (3.37) of [19] we get that in the large s limit the energy E depends on 5 n according to the 
formula 



L-X 



E = 41ns-2L^(l) + ^ 



n=2 



(B-8) 



Developing for small S n we have 

L-x r 

£ = 41ns - 2Lijj(l) + 5^ 2^ 

n=2 - 



12 



^ + 0(51 



(B.9) 



Now, for the ground state k n = L + 1 — 2n. Since we want to compare this formula with ( 12. 15ft . 
we have to consider only terms in the energy which go as In s j k and j h , with 1 < k < 5. In such 



an approximation the relevant sums appearing in (1B.9|) read as 



L-X 

E« 

n=2 



L-X 

E« 

n=2 



[2 In s-L^ (|) + (L-2)^(1)] 2 

L ^ (2) (|) 16 (^ s )2^n ~~ ¥^^(l)i6(his) 2 ^ ~~ ( 1 ) 16(ki s)2 Zlj =2 ^ 



E[fc 

n=2 



- 2- 



21ns - Lip (±) + (Z, - 2)^(1) 



si 

16 



L-X 



[21ns -L0 (i) + (L-2)#) + ...rf 



£*1 



(B.10) 



+ 



(B.ll) 
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Now, using the sums 

L-1 



£^ = ^^ + 0(L-2), J2K= { -^^ + 0(L-2f 



(B.12) 



n=2 

and also the fact that 



71=2 



21ns - Lip [ - ) + (L- 2)^(1) = 21ns 



21n2 + 7E ; 

1 H : h m 2 j 

in s 



(B.13) 



one gets, at the order of interest, 



L-1 
n=2 



L-1 
n=2 



7T 2 (L-2) J 



21n2+7B\ 2 



In s 



48(lns) 2 (1 

tt 4 ^)(f)(L-2)' 
3840 (Ins) 5 



1 - 



21n2 



3 + 



3(ln2) 5 



7T 



In s 



1280 + 21n2+ 7E j 



4 J ) 



(B.14) 
(B.15) 



where, in order to keep formula? compact, we did not develop in inverse powers of Ins the brackets 
in the denominators. With the help of the explicit values 



(2) 



-14C(3), ^ 



(i) 



-744C(5) 



(B.16) 



we eventually plug all the results in flB.9j) . getting the expression 



En 



4 In s — 4 In s In 2 j + + In s 



7<(3)x : 



24 (l + l^f 



In s 



7C(3)7r 2 ln2 



12 p | 21112+7^ ^3 



+ 



In s 



+ In s 



77r 2 (ln2) 2 C(3) 31tt 4 C(5) 
8 640 



^ Ins 



) 



+ 



49tt 4 C 2 (3) , 



960 



f + ... 



(B.17) 



Finally, developing the brackets in the denominators in their Ins and (Ins) contribution, we get 
formula (I2TT51) . 
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